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A POLAR PEAK OF ELECTRON CONCENTRATION 


AT 1000 KILOMETERS ALTITUDE 
L. II. Brace, H. G. Mayr and K. K. Mahajan 

ABSTRACT 

Satellite probe measurements from a circular, polar orbit at 1000 kilom- 
eters altitude have revealed a broad enhancement of ionization within the polar 
cap. The latitude of the maximum concentration closely follows the seasonal 
movement of the solar terminator across the polar cap. The peak is bounded on 
its nightside by £» solar zenith angle of about 110° and on its dayside by a broad 
midlatitude trough near 60° geomagnetic latitude. Thus one boundary of the peak, 
and the location of the maximum, is strongly influenced by solar zenith angle 
while the other boundary is geomagnetically fixed. Although the ionization and 
heating observed in the peak cleanly arises from solar radiation, the most in- 
triguing question is why a peak should be formed at all since this cannot be pro- 
duced by zenith angle variations alone. The answer appears to lie in the mech- 
anisms which produce the midlatitude trough on the dayside. To explain the 
trough we invoke the effects of neutral atmospheric winds flowing poleward on 
the dayside of the earth in response to solar heating. The observations are 
shown to agree quantitatively with solutions of the energy and particle equations 
of the plasma when a neutral wind field derived from a Jacchia model atmosphere 

is introduced and suitable allowance is made for an upward proton flux at high 
latitudes . 

We also conclude that the polar peak discussed here and the one observed 
by Nishida and interpreted as a neutral point phenomenon may actually be the 
same feature. Nishida 's study employed only data acquired near equinox, thus 
he could not resolve the seasonal movement we describe here. 
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A POLAR PEAK OF ELECTRON CONCENTRATION 


AT 1000 KILOMETERS ALTITUDE 


INTRODUCTION 

The high latitude topside ionosphere has been the subject of many satellite in- 
vestigations (Thomas and Sader, 1964) (Muldrew, 1965) (Thomas et al., 1966) (Nishida, 
1967) (Sato and Colin, 1969) (Jelly and Petrie, 1969) (Chan and Colin, 1969) (Brace, 1969) 
(Miller and Brace, 1969). This region exhibits a great deal of complex behavior 
in the form of peaks and troughs with considerable spatial and temporal variation. 
However several large-scale features of the high latitude ionosphere are found 
rather consistently, and these represent the background behavior of that region. 

The subject of this paper is a discussion of one of these permanent features 
called the polar peak. This feature was described by Nishida and attributed to 
solar plasma penetrating the ionosphere at the magnetospheric neutral point 
near 80° geomagnetic latitude. Our electrostatic probe measurement from Ex- 
plorer 22 do not support this explanation for the polar peak. 

THE EXPERIMENT 

The Explorer 22 satellite was launched on October 10, 1964 into a nearly 
circular orbit (1000 ± 100 km) having an 80° inclination. The cylindrical probe 
was placed aboard to provide measurements of the local electron concentration 
(N p ) and temperature (T e ). Brace and Reddy (1965) and Brace etal.,(1967) have 
described the method. The current sensitivities and voltage sweep rate employed 
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in this satellite permitted N e to be measured in th j range 1 x 10 2 — 5 x 10 4 cm -3 
and T in the range 500°— 20,000°K. The relative accuracies of the data em- 
ployed in this study are believed to be better than 10%. Comparisons with other 
probe and sounder methods (Donley et al., 1969) suggests that the absolute 
accuracy is better than 20% in N e and better than 10% for t\ . 

EXPERIMENTAL DATA 

In the Explorer 22 probe experiment, the normal mode of operation permits 
measurements to be made at 3 minute intervals along the orbit, an interval that 
permits latitudinal resolution of about 10°. The broad network of STADAN telem- 
etry stations at high northern latitudes provides nearly complete coverage of the 
polar cap. For this study data were employed from passes over England, New- 
foundland, North Central U. S., and Alaska. All of these stations can receive 
data to the highest latitude reached by the satellite (80°). 

Pole-to-Pole Plots 

To aid in surveying the T e and N (i measurements, pole-to-pole computer 
plots of each entire week of data are prepared routinely. The values are usually 
plotted against geomagnetic, invariant or dip latitude because magnetic coor- 
dinates provide a high degree of ordering of the data from a wide range of longi- 
tudes. Figure la is a typical weekly plot of N e versus geomagnetic latitude. 

The data are from the period May 12-19, 1965. 

In routine viewing of many of the weekly plots it was readily apparent that 
some features of the polar ionosphere were not as well ordered as were features 
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at lower latitudes. In particular a broad peak of N ( . within the polar cap was 
poorly resolved in geomagnetic latitude even though the concentrations in the 
peak were often as high as that found at low latitudes. A second feature, that we 
refer to here as the "midlatitude trough," alwnys marks the dayside or equator- 
ward boundary of the polar peak. Unlike the peak, this tvough is well ordered 
on a geomagnetic latitude plot and is consistently found to be centered near 60° 
geomagnetic latitude. 

The geomagnetic latitude of the peak itself was found to differ by as much as 
25’ between passes taken on the same day but at different longitudes, a fact that 
demonstrates the poor correlation of the peak with magnetically controlled 
processes. 

To investigate the possibility that the peak was produced by solar ul- 
traviolet radiation, the data in Figure la was replotted, versus geographic lati- 
tude as shown in Figure lb. The peak is notably better resolved in this plot, 
especially the location of the maximum and the steep slope which forms its 
nightside boundary. This shows that the peak is encountered by the satellite at 
the same geographic latitude regardless of the earth meridian that happens to be 
beneath the orbit at the time of the pass. In each case, as will be shown, the 
first trace of enhancement began at a solar zenith angle of between 105 and 110°, 
with the maximum concentration occurring near a zenith angle of 85°. 

Seasonal Behavior of the Peak 

Previous examination of Explorer 22 data from different seasons (Brace, 
1969) has shown that a deep trough of N p over the pole in winter gives way to 


3 




generally high values of N p in summer. Hagg (1968) reported similar results 
from Alouette-II. 

To examine this seasonal variation of the polar peak in detail, we selected 
one-week periods in winter, spring and summer and repeated this for two con- 
secutive years. In each of these periods, the satellites orbital plane was within 
3 hours of the noon-midnight plane. Thus the satellite encountered the solar 
terminator high in the polar region and crossed it nearly perpendicularly. This 
permits any solar zenith angle effects in the ionosphere to be encountered 
abruptly. Figures 2 ind 3 show the seasonal movement of the polar peak in these 
two years, 1965 and 1966. 

The initial rise in N occurs near the solar zenith angle of 105° (small 
arrow) in all six periods. The geographic latitude at this point varies widely 
with season. In summer, the 105° point occurs well into the "nightside" of the 
polar cap, and conversely it moves slightly into the "dayside” of the polar cap 
in winter. 

The weeks selected for this study were not unique. A broad survey of the 
data plots from weeks adjacent to those shown here reveals the polar peak at 
its e %pet id location. It must 1x3 realized, however, that precession of the orbit 
limits our seasonal resolution of solar terminator effects within the polar cap. 
For example, in those periods when the orbital plane tends to parallel the 
terminator, the satellite may not cross the terminator at all within the polar 
cap (When the terminator is crossed at lower latitudes w'e observe sunrise or 
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sunset effects in T e and N o that are similar in some respects to the polar peak 
enhancements but are complicated by heat and particle exchange with the over- 
lying protonosphere.) 

The Midlatitude Trough 

As noted earlier, a key feature inherent in all of the observations of the 
polar peak is the existence of the midlatitude trough of N p on its dayside. In- 
deed it is this trough that lakes the polar enhancement appear as a peak. From 
geomagnetic plots such as shown in Figure la, one finds that the location of the 
trough does not follow the seasonal movement of the peak, but remains near 60° 
geomagnetic latitude. This magnetic influence is also evident in the higher de- 
gree of ordering of the trough within the data from each week when plotted 
versus geomagnetic latitude. 

Because of the seasonal movement of the peak, its location relative to the 
midlatitude trough varies greatly. In polar summer, when the peak is furthest 
from the trough, they are separated by about 40 to 60° depending upon the longi- 
tude. Conversely, in polar winter, the peak moves into the trough region and its 
amplitude is greatly depressed. The peak may be entirely absent at European 
longitudes where the trough occurs at the highest geographic latitudes. In effect 
the mechanism that forms the trough appears to entirely suppress the peak. 

From the behavior described above it is clear that the key to the polar peak 
behavior lies in the processes that form the midlatitude trough, a topic which 
we will return to later. 
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The Electron Temperature 


If the ionization in the i>olar peak is produced by solar ultraviolet radiation, 
as seems clear from the solar zenith angle dependence just described, we can 
also expect increases in T p near the terminator. Figure 4 compares the T p and 
N ( behavior at this boundary during a 3-day period in the winter of 1965. The 
measurements trom 6 Alaska passes are plotted versus solar zenith angle. 
Typical values of N r within the winter polar trough at the right are about 500 
electrons/cc, concentrations too low to permit precise T p measuremeits with 
the current sensitivity employed on Fxplorer 22. On most passes the initial 
rise in N is found between 105° and 110° zenith angle. T p is already rising by 
the time it becomes measurable. On one pass in this series the zenith angle 
effect appears to be masked by a local enhancement of T c and N e . These occur 
si>oradically, and perhaps represent particle precipitation events (Donley, 1967). 
In any case, it is clear from Figure 4 that the large gradient in N p across the 
terminator is accompanied by corresponding changes in T e . It is significant that 
T ( . and do not exhibit a similar correlation on the dayside of the peak where 
T e continues to rise throughout the midlatitude trough of N p 

DISCUSSION 

Figure 5 is a sketch ol the relationship between the polar peak, the mid- 
latitude trough and tho solar terminator for a day in northern summer. We have 
shown no peak in the southern polar region (winter pole) assuming that its be- 
havior is similar to that of the northern pole. The limited telemetry coverage 
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of that region permits us to observe only portions of the southern polar peak 
that are equatorward of 70 ° south latitude. The geomagnetic field is represented 
by dashed lines for later discussion. 

The commencement of the N e enhancement at solar zenith angles of about 
105° is consistent with the theoretical work of Thomas (1966), although he also 
predicts increases in N e at lower altitudes that are not observed. Brinton et al. 
(1969a) reported sunrise and sunset effects beginning at similar zenith angles 
from Explorer 32 measurements at lower latitudes. Similarly the seasonal move- 
ment of the peak as it follows the terminator across the polar cap leaves little 
reason to doubt that solar illumination is intimately involved in the formation of 
the polar peak. 

It is equally clear that the trough in N p on the dayside of the peak must be 
caused by other processes, sincr. t 1 3 zenith angle continues to decrease through 
the trough region. The remainder of this paper deals with a theoretical investi- 
gation of the processes that control the midlatitude F-region and ends with a 
model which reproduces the major features of the trough and the peak. 

Qualitative Discussion 

When initially considering the origin of the midlatitude trough it is tempting 
to see if variations in temperature and composition could account for it. These 
parameters exert a strong control over the plasma scale height, H. In dif- 
fusive equilibrium H is given by 


The mean ion mass, m t , approaches 16 throughout the entire upper F-region 
where the peak is found (below 1000 km and above 60° geomagnetic latitude). T c 
exhibits a maximum in the trough region and therefore varies in the wrong 
sense to produce a trough of N c . Thus if 0 + were in diffusive equilibrium in the 
trough region, there would not be a trough unless it also existed in the lower 
F-region where ion production and loss control N p . Nishida's (1967) observa- 
tions show little evidence of the polar peak or the trough at the altitude of the 
F 2 ~max. In addition, analysis of Alouette-I sounder data from the periods of 
this study also shows little trace of the peak at that altitude. Thus we conclude 
that the trough results from some dynamic process in the atmosphere that pre- 
vents the upper F-region at midlatitude from approaching diffusive equilibrium. 
This implies the existence of a process that produces a downward force upon the 
ions that in turn reduces the electron scale heights and concentrations throughout 
the upper F-region. The observations suggest that this force is greatest at mid- 
latitudes and exerts a strong geomagnetic control. 

Kohl and King (1967) have demonstrated that the thermally driven global 
wind system can have important affects upon the global morphology of the ionos- 
phere. Volland and Mayr, (1968) have developed a model for such a wind system 
baseci on Jacchia’s model of the neutral atmosphere. Brinton et al. (1969b) have 
invoked this wand model to explain the strong hemispherical asymmetry they 
observed in the 0 + - H + transition altitude at American longitudes. Similarly 
Brace et al. (1969) showed that daytime rocket borne measurements of T t . and 
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N e within the trough region could be explained by similar wind velocities. Thus 
it appears that the midlatitude trough can be understood qualitatively in terms 
of winds in the neutral atmosphere. 

Neutral wind models show that the wind flows poleward on the sunlit side 
of the Earth. Since the ions resist movement perpendicular to the geomagnetic 
field on which they are effectively frozen, collisions between the ions and neu- 
trals induce a drag force that has a component down the field line. This pro- 
duces both a reduction of the plasma scale height and a lowering of the altitude 
and concentration at the F 2 -max. The combined effect is to reduce the electron 
and ion concentration at all altitudes above the F 2 -max. 

This process is an especially appealing mechanism for explaining both the 
midlatitude trough and the polar peak because it is most effective at midlatitudes 
where the dip of the field provides a maximum vertical component of ion-neutral 
drag as can be seen in Figure 5. At the higher latitudes where the peak is 
found, the field is nearly vertical and the largely horizontal winds have only 
small components along the field lines. This dip angle dependence is consistent 
with the remaining observational point, namely, the geomagnetic control of the 
trough location. 

The Effect of Winds on the Nightside 

Before discussing a more quantitative treatment of the midlatitude trough 
problem, it is interesting to consider the effect of the wind upon the nightside 
midlatitude ionosphere. Since the winds flow equatorward, one expects an 
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uplifting of the F-region, reduced recombination, and generally en- 
hanced pJasma concentrations at midlatitudes (Kohl and King, 1967). 
Increases in scale height further enhance the concentrations at higher altitudes. 
The nightside data in Figures la and lb tend to exhibit this kind of behavior be- 
tween 20° and 50° latitude where the electron concentration is higher than else- 
where on the nightside. Indeed, at 40° the nighttime values of N e at 1000 kilom- 
eters actually exceeds the daytime values, at least at solar minimum (Brace et al., 
1968) (King et al., 1968). 

It is puzzling that this nightside enhancement of N e is limited to lower mid- 
latitudes and is not evident at 60°. Instead there is a rather narrow trough at 
60°, called the "main trough" by Muldrew (1965), whose origin is not fully under- 
stood. Apparently a very effective depletion process is operating at that location, 
a process that reduces the electron concentration throughout the F-region (Thomas 
and Sader, 1964) (Sharp 1966). The wind-induced drag forces have only a minor 
effect by comparison. 

Analysis of the Trough and Peak Formation 

Following the method of Mayr et al. (1967) we have calculated the latitudinal 
variation of N c as a function of latitude, first for an ionosphere in diffusive 
equilibrium and then introducing processes which permit departures from dif- 
fusive equilibrium. The calculations represent a simultaneous solution of the 
particle and energy equations for the plasma (Mayr et al. 1969). They include 
the effect of the i^n a.'d electron temperatures and their gradients. Force terms 
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are included in the 0+ and H+ momentum equations to account for ion-neutral 
drag effects of the neutral wind and ion-ion drag effects of proton fluxes into 
and out of the protonosphere (Brace et al. 1969). The photoelectron flux (non- 
local heating) and the solar ultraviolet flux (local electron heating) were also 
adopted following Brace et al. (1969). 

To avoid arbitrary lower boundary conditions on N e , T e and T. , the calcula- 
tion was begun at 170 kilometers where local production and loss of 0 + estab- 
lished N e and local heating and cooling established T, and T t . T n , the gas 
temperature, was derived from the same Jacchla model as employed for the 
wind field (T ffi = 830°). The profile of T e measured at 1000 kilometers at equinox 
was employed as a second boundary condition on the temperature calculations. 
The effects of atmospheric attenuation of solar ultraviolet at low zenith angles 
was included assuming H e 304 to be the primary F-region ionization source. 

For this purpose the season was assumed to be equinox when a zenith angle of 
90° occurs at the pole. The geomagnetic and geographic poles were assumed to 
lie on a common meridian. The calculations were performed for the field lines 
30°, 45% 60°, 80° and 90° geomagnetic (and geographic) latitude, and smooth 
latitude profiles were derived by fitting those points. 

The Equilibrium Case 

Neglecting dynamic processes in the atmosphere, one expects a general de- 
crease in N e toward higher latitudes where solar zenith angle effects become 
important. To demonstrate this, Figure 6 shows calculated latitude profiles at 
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several altitudes for an ionosphere in equilibrium. No trough is present at mid- 
latitudes and the concentrations are about a factor of 4 greater than are observed. 
Clearly the observations are not consistent with the equilibrium calculations. 

The Effect of Winds 

To closer match the observations; that is, to produce a trough in N e at mid- 
latitudes, we introduced the poleward neutral wind field derived by Volland and 
Mayr (1968). Table I shows the wind speeds as a function of latitude. 

The effect of this wind on the latitudinal profiles is shown in light solid lines in 
Figure 7. At 200 kilometers and at the F 2 -peak the wind has little or no effect, 
as this region is dominated by local production and loss of 0 + . However, above 
the F 2 -peak a trough develops whose relative depth is consistent with the obser- 
vations at 1000 kilometers. But the location of the calculated trough is about 
10° too far south. 

Proton Fluxes 

In earlier studies of this type Mayr et al. (1967), Mayr et al. (1969) found 
evidence for proton fluxes flowing upward into the protonosphere at latitudes 
above 50 . Banks and Holzer (1969) also predict such fluxes associated with the 
polar wind. These proton fluxes deplete the H + population and thus decrease the 
plasma scale heights in the O 4 - H + transition region. By introducing fluxes 
that increase with latitude it is possible to shift the position of the trough to 
60 to match the observations. Figure 7 shows the resuiting solution in heavy 
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solid lines. The proton flux required is shown in Table I. These fluxes are 
everywhere less than the critical flux that can be maintained by production of 
protons through charge exchange when currently accepted hydrogen concentra- 
tions are employed (Meier, 1969). 

Other features of the theoretical trough, and the resulting peak at higher 
latitudes, are also similar to the observations. The peak itself is located at a 
zenith angle of about 80° , in acceptable agreement with the observed peak near 
85°. Our calculations are not carried to large enough zenith angles to evaluate 
the nights ide of the peak, but the observations exhibit acceptable agreement with 
the theoretical work of Thomas (1966) for large zenith angles. The growth in 
the amplitude of the peak with altitude is consistent with the observations of 
Nishida (1967). 

The Polar Peak of Nishida 

As noted earlier Nishida (1967) reported observations of a polar peak which 
he found in Alouette-I measurements taken near equinox 1963 and 1964. Nishida 
(1967) and Thomas et al. (1966) attributed this enhancement to the precipitation 
of solar wind particles entering the magnetosphere at the neutral point near 80° 
geomagnetic latitude on the dayside of the Earth. It may be less than coinci- 
dental that the polar peak we have discussed here is also found near 80° at 
equinox. Furthermore we observe no other stable peak at this location. Al- 
though the two sets of observations are not from the same years, it seems likely 
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that these two polar peaks are the same phenomenon. The seasonal movement 
of the peak may not have been readily resolvable in Nishida's data since it was 
taken near equinox. 

It should be reemphasized that other peaks are often observed in the polar 
ionosphere, and these may exhibit a high degree of magnetic ordering (Sato and 
Colin, 1969) (Colin et al., 1969) (Miller and Brace, 1969). However, these features 
are more sporadic in occurrence and often more localized in extent and should 
not be confused with the large scale phenomenon we have discussed here that is 
a consistent feature of the polar ionosphere. 

SUMMA RY 

In this study we have shown that a peak of N c exists in the upper F-region 
of the polar ionosphere, an enhancement that follows the seasonal changes in 
solar zenith angle. The peak is bounded on its nightside by zenith 
angles of 105 to 110° and is bounded on its dayside by a midlatitude trough at 
60° geomagnetic latitude. Thus the width of the peak varies with longitude. In 
winter the peak moves into the midlatitude trough region and is depressed in 
amplitude. 

The midlatitude trough and its associated polar peak can be reproduced 
theoretically by invoking ion-neutral drag forces produced by poleward winds 
flowing in the neutral atmosphere and upward proton fluxes that increase 
with latitude. The poleward winds on the dayside of the Earth produce reduc- 
tions in N p at midlatitudes where the dip of the field permits the ion-neutral 
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drag to be most effective. Conversely, equatorward winds on the nightside of 
the Earth tend to produce midlatitude enhancements of N e , a feature that is also 
observed. 

We suggest that the polar peak reported by Nishida and attributed to particle 
heating at the magnetospheric neutral point may simply represent the equinoctial 
location of the polar peak discussed here. 
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Table I 


Neutral Winds and Proton Fluxes 


Latitude 

30° 

o 

O 

50° 

60° 

70° 

00 

o 

o 

85° 

Horizontal Wind 
(m/s) 

51.6 

84.4 

11.6 

141.4 

163.8 

180.2 

187.8 

Wind Component 
IP to Field 
(m/s) 

-33.9 

-43.3 

-44.4 

-39.2 

-29.3 

-15.9 

-8.2 

Proton Flux 

0 

0 

1 x 10 7 

2.5 x 10 7 

3.5 x 10 7 

4 x 10 7 

4 x 10 


at 1000 km 
(cm -2 sec -1 ) 
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Figure lc-Pole to Pole plot of N e measurement from the week of 15-22 May 1966 versus geomag- 
netic latitude. The satellite was in a circular orbit at 1000 kilometers altitude. Individual meas- 
urements from a given pass are connected by straight lines, except when they cross the pole. 

Figure lb-Pole to Pole plot of the same N e data as in la but plotted versus geographic latitude. 
The peak and its nightside boundary are much better resolved on this scale. 
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Figure 2-Seasonal behavior of the polar peak in 1965. Each solid line represents an average of 
data from the week indicated. The arrows mark the point where the solar zenith angle is 105 . 
In all seasons the initial enhancement occurs at this location. 
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Figure 3-Seosonal behavior of the peak in 1966. As in Figure 2, the arrows 
are at a zenith angle of 105° 
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Figure 4-Relationship between T e and N # in the nightside boundary of the peak. 
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Figure 6-Calculations of the latitudinal variation of N e at several altitudes for an ionosphere in 
equilibrium. There is no evidence of rhe midlatitude trough that is observed at 1000 km (crosses). 
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Figure 7-Calculation of the latitudinal variation of N e when a poleward neutral wind is introduced 
(narrow solid line) and when a upward proton flux is added to the neutral wind (heavy line). 
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